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We consider star networks of chaotic oscillators, with all end-nodes connected only to the central 
hub node, under diffusive coupling, conjugate coupling and mean-field type coupling. We observe 
the existence of chimeras in the end-nodes, which are identical in terms of the coupling environment 
and dynamical equations. Namely, the symmetry of the end-nodes is broken and co-existing groups 
with different synchronization features and attractor geometries emerge. Surprisingly, such chimera 
states are very wide-spread in this network topology, and large parameter regimes of moderate 
coupling strengths evolve to chimera states from generic random initial conditions. Further, we 
verify the robustness of these chimera states in analog circuit experiments. Thus it is evident that 
star networks provide a promising class of coupled systems, in natural or human-engineered contexts, 
where chimeras are prevalent. 


I. INTRODUCTION 

Chimera states have been extensively studied over the 
last decade in natural and artificial networks of coupled 
identical complex systems, in fields ranging from physics 
and chemistry to biology and engineering [Ml EHSI]. 
At the outset, Kuramoto and his colleagues [THS] first 
noticed that in a system of non-locally coupled identi¬ 
cal phase oscillators, in a ring configuration, the system 
spontaneously broke the underlying symmetry and split 
into synchronized and desynchronized oscillator groups. 
Namely, there emerged a state where coherent and in¬ 
coherent sets of oscillators co-existed. This state was 
dubbed a chimera state [4], as it was reminiscent of 
the greek mythological creature composed of incongru¬ 
ous parts. Subsequently this fascinating phenomena has 
been observed in a variety of non-locally coupled sys¬ 
tems, such as time delayed systems m, Josephson junc¬ 
tion arrays m, electrochemical systems [18] and uni- 
hemispheric sleep in certain animals [19] . In recent years 
chimera states have also been observed experimentally 
in optical analogs of coupled map lattices m, Belousov- 
Zhabotinsky chemical oscillator systems m, two popu¬ 
lations of mechanical metronomes [ 20 ] and modified time 
delay electronic circuit systems m- 

Till now chimera states have been reported primarily in 
networks that have a regular ring topology, where oscilla¬ 
tors are coupled in a non-local [illllS] or global fashion 
[7| [23] . In this work we will show how chimera states also 
emerge in oscillator networks with a star topology. The 
star configuration is one where the network has a cen¬ 
tral hub position and all other nodes are linked to this 
node [16]. This configuration arises extensively in com¬ 
puter networks, where every node connects to a central 
computer, and the central computer act as a server and 
the peripheral devices act as clients. Further, a star-like 
structure is a primary motif in scale-free networks, which 
have been reported to arise in wide-ranging phenomena 

m- 

Here we will show the extensive existence of chimeras 


in the end-nodes of the star network, which are identi¬ 
cal in terms of the coupling environment and dynamical 
equations. We will demonstrate how the symmetry of 
the end-nodes is broken and co-existing groups with dif¬ 
ferent dynamical behaviour emerge. Interestingly we find 
that such chimera states are very wide-spread in this net¬ 
work topology, and large parameter regimes of coupling 
strengths typically yield a chimera state. We also confirm 
the existence of robust chimera states in analog circuit 
experiments. 


II. STAR NETWORKS OF CHAOTIC 
OSCILLATORS 

Here we study the dynamics of a star network of N 
identical nonlinear oscillator systems. In such networks 
there is one central hub node (labelled by site index i = 1) 
and N — 1 environmentally identical peripheral end-nodes 
connected to the central node (labelled by node index 
i = 2,... N). One can also interpret this system as a set 
of uncoupled oscillators connected to a common drive. 
The focus of this study is the dynamical patterns arising 
in the N — 1 identical end-nodes of this network. In or¬ 
der to establish the generality of our results, we consider 
three different coupling forms: (a) diffusive coupling (b) 
conjugate coupling and (c) mean-field coupling. We give 
below the general dynamical equations for the different 
coupling forms. First, we consider standard diffusive cou¬ 
pling through similar variables, given by: 

N 

~ fyit ^ ^ ^ij ( 1 ) 

i=i 

Vi = fy{Xi,yi,Zi) 

— fz{xi^ yi^ Zi) 

Here coupling matrix element for central node i = 1 
is Kij = k/2 when j 7 ^ I, and for the end-nodes 
i = 2,... Kii = /c/2 and zero otherwise. The coupling 
strength is given by k. Then we consider the conjugate 


coupling [24] given as: 
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N 

Xi = fx{xi,yi, Zi) + ^ Kij{yj - Xi) (2) 

i=i 

m = fy{Xi,yi,Zi) 

Zi = fz{xi,yi,Zi) 

Lastly, we also consider a mean-field type of coupling, 
where the dynamics of the central node is given by: 

k 

Xi = fx{xi,yi,Zi) + -{Xm-Xi) (3) 

m = fv{xi,yi,zi) 

Zi = fz{xi,yi,zi) 

where Xm = Sj =2 n mean field of the 

end-nodes. The dynamics of the end-nodes i = 2,... N 
is given by: 


k 

Xi = fx{xi,yi,Zi) + -{xi - Xi) (4) 

m = fy{Xi,yi,Zi) 

Zi = fz{xi,yi,Zi) 

For the local dynamics at the nodes, we take two pro¬ 
totypical chaotic systems that have widespread relevance 
in modelling phenomena ranging from lasers to circuits. 
First we consider the Rossler type oscillator at node i, 
given by the form: 

fx{xi,yi,Zi) = -[uji + e{xi‘^ ^ yi‘^)]yi - Zi (5) 

fy{xi, Pi, Zi) = [cji + c{xi^ -h y?)\xi + ayi 

fz{xi,yi,Zi) = b + Zi{xi -c) 

in Eqns. 1-4. For each node, uJi + e{xi‘^ + yi^) is close 
to the angular velocity of the oscillator, perturbed 
by amplitude Xi‘^ + yi^ when e 7 ^ 0. Here we take the 
parameter values to be: a = 0.15, 6 = 0.4, c = 8.5, 
cji = ci ;2 = Co’s = 0-41 and e = 0.0026 [24|, yielding a 
chaotic attractor. We also consider the Lorenz system at 
the nodes given by: 

fx{xi,yi,Zi) = (j{yi - Xi) ( 6 ) 

fy{xi, yi, Zi) = (r - Zi)xi - yi 

fz{xi,yi,Zi) = Xiyi - f3zi 

in Eqns. 1-4. With no loss of generality we consider the 
parameters of the local system to be a = 10 , r = 28 
and /3 = 8/3, yielding double-scroll attractors. We study 
both these chaotic systems, coupled in star network con¬ 
figuration, through the different coupling forms given 
above. A wide range of coupling strengths, in networks of 
size ranging from 3 to 100 oscillators is investigated. The 
principal observations of the patterns arising in these net¬ 
works, from generic random initial states, are described 
below. 


A. Dynamical Patterns for Coupled Rossler 
Oscillators 
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Figure 1. Time evolution of the x variable for oscillators in 
distinct synchronized groups in a coupled Rossler systems in a 
star network of (a) 10 conjugately coupled nodes (cf. Eqn. 2) 
for coupling strength /c = 0.22; (b) 100 conjugately coupled 
nodes (cf. Eqn. 2) for coupling strength k = 0.18; (c) 10 
diffusively coupled nodes (cf. Eqn. 1) for coupling strength 
k = 0.14; and in (d) 100 diffusively coupled nodes (cf. Eqn. 
1) for coupling strength /c = 0.14. In all these figures the 
identical end-nodes, that split into 3 groups in (a) and (b), 
and into 2 groups in (c) and (d) are marked with dotted and 
dashed lines of different colors. The one distinct central hub 
node is also shown (in solid red) for comparison. 


We find that as coupling strength increases, the end- 
nodes go from a de-synchronized state to a completely 
synchronized state, via a large coupling parameter regime 
yielding chimera states. In the representative examples 
of chimera states displayed in Eig. the 9 identical end- 
nodes of the network of 10 conjugately coupled oscillators 
(Eig ^), clearly split into 3 clusters, with two synchro¬ 
nized clusters having 4 oscillators each, and 1 oscillator 
being distinct from both these synchronized groups. Eor 
the case of 100 conjugately coupled nodes in Eig. [^, the 
99 identical end-nodes again cluster into 2 synchronized 
groups of size 49 each, and 1 oscillator is uncorrelated to 
either group. Notice that the central node settles down 
to low amplitude oscillations, while the end-nodes ex¬ 
hibit large amplitude oscillations, with each group hav¬ 
ing a different phase with respect to another. Eor the 
coupling strengths presented in the figure, regular low- 
period oscillations emerge in the end-nodes, though the 
constituent oscillators were chaotic. 

Eor the case of 10 diffusively coupled Rossler oscilla¬ 
tors, the identical end-nodes split into 2 synchronized 
clusters of sizes 5 and 4 (Eig[^), while the end-nodes of 
a star network of 100 diffusively coupled oscillators split 
into 2 synchronized clusters of size 53 and 46 oscillators 
(Eig[^). Here the central node and the end-nodes all 
exhibit large amplitude oscillations of higher periodicity. 

Fig-i shows the oscillatory patterns of the end-nodes 
for the distinct synchronized groups that emerge from 
generic random initial states. Eor instance, it is evi- 
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dent from Fig. (a-d) that sub-sets of the end-nodes 
display very different attractor geometries, though they 
have identical dynamical equations. So from Figs. 
and it is clearly evident that chimera states emerge 
in the end-nodes of the star network. Further, Fig. 
shows the state of synchronization of the different end- 
nodes i = 2 ,... A/" at some representative instant of time, 
demonstrating the co-existence of synchronized and de¬ 
synchronized groups among the identical N—1 peripheral 
nodes in the star network. Note that there is no space 
ordering of the node index i of the end-nodes. So the 
(de)synchronized nodes in a cluster are not “contiguous”, 
as is usual in regular lattice topologies. 
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Figure 3. A matrix displaying the state of synchronization 
of nodes i and j in a star network of conjugately coupled 
Rossler systems {i,j = 2,...N). The blue color indicates 
that the nodes are synchronized and the green that they are 
desynchronized. Here coupling strength k = 0.24 and system 
size N = 100. The presence of a synchronized group of nodes, 
along-side a desynchronized set, can be clearly seen. 



B. Dynamical Patterns for Coupled Lorenz systems 




Figure 2. Phase portraits of conjugately coupled Rossler sys¬ 
tems in a star network with (a) 10 nodes and coupling strength 
k = 0.12, (b) 100 nodes and coupling strength k = 0.22. 
Phase portraits of diffusively coupled Rossler systems in a star 
network with (c) 10 nodes and coupling strength k = 0.14, (d) 
100 nodes and coupling strength k = 0.14. In all these fig¬ 
ures the distinct end-node clusters are marked in blue and 
cyan. The central hub node is also shown (in solid red) for 
comparison. 


Here again we find that as coupling strength increases, 
the end-nodes go from a de-synchronized state to a com¬ 
pletely synchronized state, via a large coupling parameter 
regime yielding chimera states. We display some repre¬ 
sentative patterns from the chimera states in Fig. It 
is clearly evident from these that the identical end-nodes 
split into different dynamical groups, thereby breaking 
symmetry. Some of these groups consists of synchronized 
nodes and some are clusters of de-synchronized elements, 
as seen from Fig. Further, it is also evident from Fig. 
I^that in addition to different synchronization properties, 
the groups also yield different attractor geometries. 



Figure 4. Temporal patterns (a-b) and their corresponding 
phase portraits (c-d) for a star network of coupled Lorenz 
system with (a) 10 conjugately coupled nodes and coupling 
strength k = 2.9, yielding 2 synchronized clusters of size 6 
and 3; (b) 10 diffusively coupled nodes for coupling strength 
k = 13, yielding one synchronized cluster of size 7, and a 
desynchronized cluster of size 2. 
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Figure 5. State Xi of the end-nodes of a star network of con- 
jugately coupled Lorenz systems (z = 2,... N)^ at an instant 
of time. Here coupling strength k = 1.98 and system size 
N = 100. The presence of a synchronized group of nodes, 
along with a desynchronized set, can be clearly seen. 



Figure 7. Synchronization error of the end-nodes of a star 
network of diffusively coupled Lorenz systems as a function 
of time (namely the standard deviation of i = 2,... A/", 
at an instant of time). Here coupling strength k = 5.12 and 
system size N = 100. It is clearly evident from the oscillating 
synchronization error that the end-nodes move in and out of 
synchronization. 


Further we find that the incoherent state may be 
desynchronized at the same level (stable chimera) or yield 
an oscillating incoherent group which goes in and out of 
synchronization, namely a breathing chimera [27]. Such 
a breathing chimera state is displayed in Figs. IZl The 
occurence of breathing chimera states is more common in 
the coupled Lorenz system than in coupled Rossler sys¬ 
tems. In fact breathing chimeras were also observed in 
Lorenz systems coupled in a ring configuration in earlier 
studies |9|. 
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Figure 6. Temporal patterns of the end-nodes of a star 
network of diffusively coupled Lorenz systems displaying a 
breathing chimera state. Here coupling strength /c = 5.0 and 
system size N — 100. 


C. Prevalence of chimera states 


In order to quantify the probability of obtaining 
chimera states from random initial states we calculate the 
fraction of initial conditions leading to co-existing syn¬ 
chronized and desynchronized states in the end-nodes, in 
a large sample of random initial states. This provides an 
estimate of the basin of attraction of the chimera state, 
and indicates the prevalence of chimeras in this system. 
So this measure is important, as it allows us to gauge the 
chance of observing chimeras without fixing special initial 
states. 

Figs. and display this quantity for star networks 
of Rossler and Lorenz systems. It is clearly evident from 
these figures that there exists extensive regimes of cou¬ 
pling parameter space where the probability of obtaining 
a chimera state is close to one. This quantitively estab¬ 
lishes the prevalence of chimeras in the end-nodes of non¬ 
linear oscillators coupled in star configurations. Also no¬ 
tice that larger networks yield larger basins of attraction 
for the chimera state. Further, the figures show that con¬ 
jugate coupling yields larger parameter bands with high 
prevalence of chimera states. 
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Figure 8. Probability of obtaining chimera states (red), 
synchronized clusters (magenta), fully synchronized states 
(green), and completely de-synchronized states (blue) in star 
networks of coupled Rossler systems, for the following cases: 
10 nodes under (a) conjugate coupling and (b) diffusive cou¬ 
pling; 100 nodes under (c) conjugate coupling and (d) diffusive 
coupling. 



\ 

\ 

K 


\ 

f 

/ 

¥ . 
f . 


0 0.5 1 1.5 2 2.5 3 3.5 0 0.5 1 1.5 2 2.5 3 3.5 

k k 


Figure 9. Probability of obtaining chimera states (red), 
synchronized clusters (magenta), fully synchronized states 
(green), and completely de-synchronized states (blue) in star 
networks of coupled Lorenz systems, of 100 nodes, under (a) 
conjugate coupling and (b) regular diffusive coupling. 


III. EXPERIMENTAL VERIFICATION OF 
CHIMERA STATES 


Now we establish the robustness of these chimera states 
in experimental situations by demonstrating the occur¬ 
rence of chimera states in star networks of coupled non¬ 
linear oscillators, evolving from generic initial states. 
Specifically, we consider a circuit implementation of a 
chaotic Rossler-type oscillator at the nodes, represented 
by the equation [25] : 


(Px 



dx 

dt 


±(N-i) 


( 7 ) 


This equation illustrates a jerk type chaotic system [25] . 
and an analog simulation circuit of this equation can 
be carried out with standard operational amplifiers and 
diodes. The details of a straightforward circuit imple¬ 
mentation of Eqn. can be found in Ref. [25] . We then 
go on to set up 4 diffusively coupled oscillators, with pa¬ 
rameter A = 0.58 in Eqn. such that the oscillators 
individually exhibit chaotic dynamics. Specifically, we 
experimentally study the star network given schemati¬ 
cally in Eig. where the central node evolves as: 

k 

Xi =yi^ -{x2-\-xs-\-X4- 3xi) (8) 

yi = 

ii = -yi± {\xi \ - 1) 

The evolution of the identical end-nodes is given by: 
k 

Xi = yi^ - {xi -Xi) (9) 

yi — 

Zi = -Azi -yi± {\xi \ - 1) 
where i = 2, 3,4 and k is the coupling co-efficient. 
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Figure 10. Schematic of the star network realized through 
analog simulation circuits. Here block 1 represents the central 
node and blocks 2,3, and 4 represent end-nodes. 


Fig. [TT] depicts the electronic analog circuit to imple¬ 
ment the central node or end-nodes of Eqns. [9][TQj If the 
circuit of Fig. acts as a central node, then the out¬ 
put voltages from op-amp U1 and U2 correspond to —xi 
and xi of Eqns. 8-9. If we use the circuit of Eig. [TT]for 
the end-nodes, then they generate the signals V2 = X 2 , 
V3 = xs and V4 = X 4 at output of 1/2. Eor the central 
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node circuit, the input voltage Vc is the coupling voltage 
signal generated from the circuit of Fig. In this case, 
the Vc signal corresponds to {k/ 2 ){x 2 + X 3 + X 4 — 3xi). 
Fig. depicts the circuit used to implement the cou¬ 
pling between central node to the end-nodes. In this case, 
Vc signal corresponds to {k/2){xi — Xi). Representative 



Figure 11. Circuit implementation of Eqns. 8-9 using op- 
amps {AD712 or fiA74:l). The external voltage supply (FI) 
is +1.0F. The capacitors are lOnF. Resistors R4 and R5 = 
l/cQ, the variable resistor is 10/cQ, and the rest are lOOkQ. 
The diode is IN4148. Vc is the input coupling signal. 
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Figure 12. Circuit used to implement coupling between the 
end-nodes and the central node. All op-amps are AD712 or 
liA741. Resistor R1 = 100/cQ, and the rest are SOOkQ. Vc 
is the output coupling voltage, and FI = —xi is signal gen¬ 
erated by the central node. Signals F2, F3 and F4 corre¬ 
sponding to X 2 , X 3 and X 4 , are generated by using three more 
circuit copies of Fig.ll. 

circuit simulation results are displayed in Fig. where 
phase-portraits in the Xi — yi plane are displayed for dif¬ 
ferent coupling strengths k. One clearly notices that for 
low coupling strength (e.g. k = 0.1 in FigpA^) the end- 
nodes show completely unsynchronized oscillations. For 
large coupling strength (e.g. k = 2.0 in Figjl^), as antic¬ 
ipated, the end-nodes exhibit complete synchronization. 
However, for moderate coupling strengths (e.g. k = 1.0 in 


R R kR 



Figure 13. Circuit used to implement coupling between the 
central node and the end-nodes. All op-amps are AD712 or 
IjlA741 and the resistors are lOO/cfl. Fc is the output coupling 
signal, Vj = xi from the central node and Vi = —xi from the 
end-nodes. 

Figjl^) the 3 identical end-nodes split into two groups, 
where two of them are synchronized and one is not, thus 
exhibiting a chimera-like state. The time series of this 
state is shown in Fig. [^to further illustrate the broken 
symmetry of the three identical end-nodes in the star net¬ 
work. Note that we have no control over the initial state 
in the experiment, and these states evolve from generic 
random initial conditions. 
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Figure 15. Times series of the 3 end-nodes (a-c) of 4 diffusively 
coupled chaotic Rossler-type circuits, represented by Eqn 
with coupling strength k — 1.0. Clearly 2 nodes (a,b) are 
synchronized, while one (c) is distinct from this group. 

Further, in order to check the generality of the re- 
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Figure 14. Phase portrait of attractors in the xi — yi plane, 
generated from chaotic Rossler-type oscillator circuits repre¬ 
sented by Eqn. that are diffusively coupled in a star net¬ 
work, for coupling strengths: (a) A; = 0.1 yielding an unsyn¬ 
chronized state, (b) k = 1.0 yielding a chimera-like state and 
(c) k — 2.0 yielding a synchronized state. In all these fig¬ 
ures the oscillations of the end-nodes in the star network are 
marked in different colours. The central hub node oscillator 
dynamics is not shown. 


suits, we also investigate the mean-field type of cou¬ 
pling given by Eqns. 3-4. Fig. displays represen¬ 
tative phase-portraits in the Xi — yi plane for different 
coupling strengths k. Again one finds that for low cou¬ 
pling strengths (e.g. k = 0.02 in Figjl^) the end-nodes 
are completely unsynchronized, while for high coupling 
strengths (e.g. k = 2.0 in Fig|T^) they are completely 
synchronized. However, in a large window of moder¬ 
ate coupling strengths (e.g. A: = 1.0 in Figjl^) the 3 
identical end-nodes split into two groups, where two of 
them are synchronized and one is not, thus exhibiting 
a chimera-like state. Also note the different geometries 
of the dynamical state in the two groups. Lastly, we 
estimate the probability of obtaining the chimera state 
in the star network with mean-field coupling by find¬ 
ing, through numerical simulations, the fraction of ini¬ 
tial states that evolve to chimera states. The results are 


displayed in Fig. and it is clear that this form of cou¬ 
pling yields a large parameter regime where the typical 
initial state gives rise to a chimera state in the end-nodes. 



Figure 16. Phase portrait of attractors in the Xi — yi plane, 
generated from chaotic Rossler-type oscillator circuits repre¬ 
sented by Eqn. coupled via mean-field in a star network 
(cf. Eqns. 3-4), for coupling strengths: (a) /c = 0.1 yielding 
an unsynchronized state, (b) A: = 1.0 yielding a chimera-like 
state and (c) k = 2.0 yielding a synchronized state. In all 
these figures the oscillations of the identical end-nodes in the 
star network are marked in different colours. The central hub 
node oscillator dynamics is not shown. 






























Figure 17. Probability of obtaining chimera states (red), 
synchronized clusters (magenta), fully synchronized states 
(green), and completely de-synchronized states (blue) in star 
networks of coupled Rossler systems with mean-field type cou¬ 
pling (cf. Eqns. 3-4) for a network of 100 nodes. 


IV. CONCLUSIONS 

In summary, we have investigated star networks of dif¬ 
fusively and conjugately coupled nonlinear oscillators, 
with all end-nodes connected only to the central hub 
node. Though the end-nodes are identical in terms of 


the coupling environment and dynamical equations, they 
yielded chimera states. Namely, the symmetry of the 
end-nodes was broken and co-existing groups with dif¬ 
ferent synchronization features and attractor geometries 
emerged. We estimated the basin of attraction of chimera 
states by evaluating the fraction of initial states that 
evolve to a chimera state, in a large sample of random 
initial conditions. This measure showed that in exten¬ 
sive regimes of coupling parameter space the probability 
of obtaining a chimera state is close to one. Further, we 
established the robustness of these chimera states in ana¬ 
log circuit experiments. The experimental verifications 
incorporated both diffusive coupling and mean-field type 
coupling for the central node. Thus it is clearly evident 
from our numerical and experimental investigations that 
large parameter regimes of moderate coupling strengths 
yield chimera states from generic random initial condi¬ 
tions in this network topology. So star networks provide a 
promising class of coupled systems, in natural or human- 
engineered contexts, where chimeras are pervasive. 
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